Moment Approach for Determining the Orbital Elements 
of an Astrometric Binary with Low Signal-to-noise Ratio 



43 



Hirofumi Iwama, Hideki Asada, Kei Yamada 
Faculty of Science and Technology, Hirosaki University, Hirosaki, Aomori 036-8561 

asada@phys. hirosaki-u. ac.jp 

(Received ; accepted ) 
^ . Abstract 

A moment approach for orbit determinations of an astrometric binary with low 
signal-to-noise ratio from astrometric observations alone is proposed, especially aim- 
ing at a close binary system with a short orbital period such as Cyg-Xl and also at 
a star wobbled by planets. As an exact solution to the nonlinearly coupled equation 
system, the orbital elements are written in terms of the second and third moments 
of projected positions that are measured by astrometry. This may give a possible 
estimation of the true orbit. 
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1. Introduction 

Space astrometry missions such as Gaia and JASMINE are expected to reach a few micro 
arcseconds (Mignard 2004 ; Perryman 2004 ; Gouda et al. 2007 ). Moreover, high-accuracy 
VLBI is also available. 

For visual binaries, formulations for orbit determinations have been well developed since 
the nineteenth century (Thiele 1883 ; Binnendijk 1960 ; Aitken 1964 ; Danby 1988 ; Roy 1988 
). At present, numerical methods are successfully used (Eichhorn and Xu 1990 ; Catovic and 
Olevic 1992 ; Olevic and Cvetkovic 2004 ). Furthermore, an analytic solution for an astrometric 
binary, where one object can be observed and the other such as black holes and neutron stars 
is unseen, has been found (Asada et al. 2004 ; Asada 2008 ). The solution requires that 
sufficiently accurate measurements of the position of a star (or the photo-center of a binary) 
are done more than four times for one orbital period of the binary system. 

For a close binary system with a short orbital period, we have a relatively large un- 
certainty in the position determination. For instance, the orbital period of Cyg-Xl is nearly 
6 days, which are extremely shorter than that of normal binary stars, say a few months and 



1 



several years. Because of such an extreme condition, it is interesting to seek another method in 
addition to the standard one. Moreover, stars with planets also are another interesting target. 

What can we do for orbit determination from position measurements with low signal-to- 
noise (SN) ratio? It is expected that the position of the object is measured many times. The 
dense region of the observed points is corresponding to the neighborhood of the apastron of the 
Kepler orbit, because the motion of the source star is slower according to the Kepler's second 
law. On the other hand, a region of fewer points is including the periastron, around which 
the source star moves faster. Therefore, a statistical analysis including the variance of the 
measured positions and their correlation will bring the information about the orbital elements 
of the binary system. Gedanken experiments suggest that the second moments are useful for 
exploring the shape of the orbit but they are not sufficient for the full orbit determination. At 
least the third moments seem to be needed. See Figs. 1 and 2. 

Therefore, the main purpose of this paper is to propose a method for orbit determination 
of an astrometric binary with low SN ratio by using the second and third moments. We shall 
provide also an exact solution for the coupled equations. As a result, the orbital elements of 
the binary are written in terms of the second and third moments. 

This paper is organized as follows. We present a formulation and the solution in § 2. In 
§ 3, numerical tests are also done to see how reliable the analytic result is for practical cases 
taking account of observational noises. § 4 is devoted to Conclusion. Throughout this paper, 
the spatial coordinates are the angular positions normalized by the distance to the celestial 
object. 

2. Basic Formulation 

2.1. Motion in the orbital plane 

Before projecting onto the celestial sphere, we consider a Kepler orbit, where the semi- 
major axis and semiminor one are denoted as clk and bx- We choose the (X,Y) coordinates 
on the orbital plane such that the X-axis is along the semimajor axis of the elliptic orbit, the 
Y-axis is along the semiminor one and the origin of the coordinates is chosen as the center of 
the ellipse. The orbit is expressed as 



For an object in the Kepler motion, the time t is related with the position u through 
the so-called Kepler equation 



X = ok cos u 



Y = bx sin u 



(1) 
(2) 



where u denotes the eccentric anomaly. 

By introducing the eccentricity ck, the semiminor axis becomes 




(3) 
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t = t Q + —{u- e K smu), (4) 

where to denotes the time of periastron passage. This equation represents transcendental equa- 
tions in the sense that it cannot be solved analytically without any approximation. In other 
words, u cannot be expressed by using elementary functions of t. This makes the inverse 
problem in astrometry difficult. 

2.2. Projection with respect to the line of sight 

We assume a plane perpendicular to the line of sight, where this plane is safely considered 
a small part of the celestial sphere because the angular position shift of extra-solar objects is 
sufficiently small. 

A reference frame, which is used in astrometric observations, is described by the (x,y) 
coordinates. As usual, the (X,Y) coordinates on the orbital plane and the (x,y) coordinates 
on the reference plane are related through 

x = (X — aKek)(cosu cosQ — sin uj sin VL cos i) 



— y(sinu;cosf2 + cos uj sin Q cos i) 

= xo + acosw + j3sinu, (5) 
y = (X — a K^k) (cos u sin Q + sin uj cos Q cos i) 
— Y(sinu;sinf2 — coso>cosf2cosi) 

— Vo + 7 cos u + <5sinw, (6) 
where x , yo, a, P, 7, S are defined as 

xo = —aKexicosucosfl — sin uj sin Q cos i), (7) 

yo = —aKexicosuismfl + sinwcosficosz), (8) 

a = clk (cos uj cos Q — sin uj sin Q cos i), (9) 

f3 = — &ft-(sina>cosf2 + cos uj sin Q cos i), (10) 

7 = a^(cosa;sinf2 + sin uj cos Q cos i), (11) 

5 = — 6/<(sina;sinf2 — cos uj cos Q cos i). (12) 



Here, f2, uj and i denote the longitude of ascending node, the argument of periastron and the 
inclination angle, respectively. See Fig. 3. 

The ascending node and the descending one cannot be distinguished by astrometric 
observations alone. Therefore, both of ± for i are possible. Furthermore, this paper focuses 
on the moments, so that the clockwise and ant i- clockwise motions cannot be distinguished. 
Two pairs of (Q,uj) are possible. Nevertheless, the shape of the orbit is uniquely determined as 
shown below. 
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2.3. Moment Formalism 

Let us assume frequent observations of the angular position. Namely, we consider a large 
number of observed points. For such a case, the statistical average expressed as a summation 
is taken as the temporal average in an integral form as 



where < > denotes the mean and T b s denotes the total time duration of the observations. 

In this paper, we focus on the periodic motion, so that the above expression becomes 
the integration over one orbital period. We thus obtain 

1 fto+P 



where we used Eq. (4) and dt = (1 — eKCOsu)du. 

Let us consider statistical moments. Figure 1 suggests that the moments M xx and M yy 
are useful for distinguishing two different orbits. Next, we consider the periastron for the orbit 
denoted by the closed solid curve. For simplicity, we assume that it is located at positive x (in 
the right hand side of the ellipse). For this case, the object moves faster around the periastron 
and slower around the apastron. Hence, the dots in the figure schematically show asymmetry 
in the number of observed points. In order to distinguish such an asymmetry, the third moment 
such as M xxx seems useful because of the odd parity. Figure 2 suggests that M xy is needed via 
a Gedanken experiment. The second moments are defined as 




(13) 




(14) 



M xx = < (x- < x >) 2 > 



< 2 , o2\ 1 2 2 

= -{a +(3 )--e K a , 



(15) 



M yy = <(y-<y >) 2 > 



(16) 



= < (x- < x >){y- <y>)> 
= ^1 + P6) --e 2 K aj, 



(17) 



The third moments are defined as 



M, 



xxx — 



yyy — 



< (x— < x >) 3 > 
^e K a(a 2 + (3 2 )-^e 3 K a 3 , 

< (y-<y>) 3 > 
i^7(7 2 + 5 2 )-^e| 7 3 , 



(18) 



(19) 



xxy — 



< (x- < x >) 2 (y- <y>)> 
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= ^(3a 2 7 + /? 2 7 + 2a/M)-^a 2 7 , (20) 
8 4 

M xyy = < (x- < x >){y- <y>f> 

= le K (3aj 2 + aS 2 + 2/3 7 <5) - ]e 3 K a^ 2 . (21) 
8 4 

The moments M xx , ■ ■ ■ ,M xyy are actually observables. 

The vanishing of all the third moments leads to ex = 0. In the following, let us consider 
ex 7^ cases. 

In the above second and third moments, the specific combinations among a, {3, 7, 5 
frequently appear. Therefore, it is convenient to define new variables as 

h = o? + p, (22) 
h = l 2 + 5\ (23) 
I 3 = a-y + p6. (24) 

By using the moments expressed by Eqs. (15)-(17), the variables I\, I2, I3 are rewritten as 

h = \e 2 K a 2 + 2M xx , (25) 

h = \eW + 2M yy) (26) 

h = ^e 2 K a 1 + 2M xy . (27) 

Eq. (25) is substituted into Eq. (18) to obtain a cubic equation for e^a as 

(e K a) 3 -12M xx (e K a) + 16M xxx = 0. (28) 

This cubic equation gives three roots as e^a by using Cardano's formula. For saving the space, 
we do not write down the formula. In the similar manner, we obtain from Eqs. (26) and (19) 

(e Kl f - 12M yy (e Kl ) + 16M yyy = 0. (29) 

This provides ex7, where its multiplicity is three. 

At this point, we know both e^a and ex7, which can be substituted into Eqs. (25)-(27). 
We thus obtain the value of the variables I±, I2 and I3. 

In general, we get multi values for e K a and 7. However, they are true solutions for 
Eqs. (28) and (29), whereas not all of them satisfy the observed third moments. They have to 
satisfy the remaining set of the third moments as 

M xxy = -{e K j)Ii + \{e K a)h - \{eKOt) 2 {e K l)i (30) 
8 4 4 

M xyy = ^{e K a)I 2 + ^(e^7)/ 3 - ^(e^a)(e^7) 2 , (31) 

which can be used to pick up the correct e^a and exj from multiple candidate values. 
The definition of I\ by Eq. (22) is rewritten as 
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which tells us a as a function of f3 

a = ±yJ-/3 2 + h. (33) 

It is obvious that the sign of the right hand side of Eq. (33) must be the same as that 
of the left hand side, namely a. For ck 7^ 0, the sign of a is the same as that of e^a that has 
been obtained above. This is expressed by 

sgn(a) = sgn(e^a), (34) 

where sgn denotes the sign. Therefore, the sign of the right hand side of Eq. (33) is obtained 
uniquely as sgn(exa). Eq. (33) thus becomes 



a = sgn(e K a) \J -f3 2 + I x . (35) 

What is the difference between Eqs. (33) and (35)? Eq. (33) means two different equations 
because of ± in the right hand side, whereas Eq. (35) is a single one. 
For the later convenience, we define V as 

r = l, (36) 

a 

which is obtained from known quantities e^a and exl as T = (e^7)(e j R-a!)~ 1 . 

We substitute Eq. (36) into the definition of ^3 by Eq. (24) in order to delete 7. We 

obtain 

2 ^3 — /5<5 

a = — - — . (37) 
This is substituted into the definition of I\ by Eq. (22). We obtain 

«=^. (38) 

This is a function of only f3. 

Eqs. (36) and (38) are substituted into 7 and S in Eq. (23). After rather lengthy 
calculations, we get 

2 _ ( J 3 - hT) 2 

p ~ hv 2 +i 2 -2ir { } 

where we used Eq. (22) for a 2 . Interestingly, the right hand side of this equation consists 
of only the known quantities I x , I2, I3 and T. Therefore, we obtain the value of j3 2 , which 
determines a through Eq. (35). 
Eq. (39) is solved for (3 as 

(3 = ± - h ~ hY ■ (40) 

Unfortunately, we do not know ex/3 contrary to e^a. Therefore, the sign of the right hand side 
of Eq. (40) is not uniquely determined. The multiplicity of (3 is two and hence that of 5 is also 
two according to Eq. (38). 
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Up to this point, we know the value of a and exot from observed quantities. We thus find 
separately ex as ex = (e^c^aT 1 . It is crucial in the following procedure that the eccentricity 
ex is determined at this step. 

It is very inconvenient that j3 and 5 are proportional to the semiminor axis bx in their 
definition. We know ex, so that bx can be expressed by ax- Hence, we define renormalized 
quantities as 



5 = 



\-e\ 

— (sin w cos + cos a; sin cos i), (41) 
5 



1-e? 



K 

= — aft-(sinu;sinf2 — cos a; cos cos i), (42) 

where we used Eqs. (10) and (12). We know the values of (3, 5 and ex- Therefore, one can 
estimate $ and 5. 

One can construct from four variables a, f3, 7, 5, some quantities that are dependent on 
the inclination angle i but not on any other angles u nor Q. One example is 

= <4(1 + cos 2 i). (43) 

Another is 

D = a5 - /3 7 

= a 2 K cosi. (44) 

These relations can be verified by direct calculations. Note that D must be positive because 
cosz > 0. This positivity chooses one pair of ((3,5) and reject the other pair. Only the pair of 
(f3,5) is thus obtained. We know a, 7, (3, 5, so that C and D can be estimated. 
By deleting a\ from Eqs. (43) and (44), we obtain 
C 

cos 2 z — — cosi+ 1 = 0. (45) 

This is a quadratic equation for cosz. Apparently, two cases of cosi are possible. However, this 
is not the case. By using Newton's identities for Eq. (45), we obtain the identity as 

(cosii) x (cos^) = 1, (46) 

where cosii and cosz2 denote the two roots for Eq. (45). The inequality as 

I coszj < 1, (47) 

thus leads to the unique value of cosi, because the other exceeds the unity as its absolute value. 

The single value of cosz provides positive i and negative one. They are corresponding 
to the ascending node and the descending node, respectively. They cannot be distinguished by 
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astrometric observations alone. In order to distinguish them, radial velocity measurements for 
instance are needed. 

We get the value of the inclination angle. From Eq. (44), therefore, the semimajor axis 
is obtained as 



a K = J—. (48) 
V cosz 

This suggests that the value of ax is uniquely determined but not doubly. 

In order to determine u and Q, let us consider other combinations among a, 7, j3, 6. 
Direct calculations lead to 

a 2 + (3 2 = a^(cos 2 fi + sin 2 ficos 2 z), (49) 
7 2 + (5 2 = a|-(sin 2 fi + cos 2 ficos 2 z)- (50) 

The ratio of them is denoted as 

7 2 + ^ 2 
r 1 = — 

a 2 + (3 2 

sin 2 Q + cos 2 Q cos 2 % , , 

— (51) 

cos 2 f2 + sin 2 CI cos 2 % 

This is solved for Q as 

9^ r 1 — cos 2 i 

t a n 2 fi = -i -, 52 

1 — n cos^ z 

which gives the values of Q because we have already determined i and r%. As mentioned already, 
both of ±f2 are allowed. 

Next, we consider different combinations as 

a 2 + 7 2 = a 2 K (cos 2 u + sin 2 a; cos 2 i) , (53) 
/3 2 + 5 2 = a 2 K (sin 2 tu + cos 2 tocos 2 i), (54) 

which can be verified by direct calculations. The ratio of them is denoted as 

P 2 + ~5 2 
a z + 7^ 

sin 2 oj + cos 2 to cos 2 i 

= — 9 — ; — ~ ( 55 ) 

cos z w + sm to cos z 1 
This is solved for to as 

2 r 2 -cos 2 z 

tan u=- —, (56) 

1 — r 2 cos^ z 

which tells us the values of u because we have already determined % and r 2 . As mentioned 
already, both of ±tu are possible. 

At most four values of Q are possible. Similarly, the maximum multiplicity of u is four. 
In total, sixteen sets of appear to exist. The multiplicity of (u,Q) is reduced, because 

they must satisfy the definition of a, (3, 7, 5 with the uniquely determined ax, £k and cosi. 
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In particular, a, f3, 7, 5 include different combinations of sin and cos. Basically, the sign of sin 
and cos has four types (+,+), (+,—), (—,+), (—,—). Hence, the apparent multiplicity sixteen 
is divided by four. In addition, it is obvious from Eqs. (9)-(12) that the sign of Q depends 
on that of u. The multiplicity is thus divided also by two. As a consequence, we get only the 
two pairs of One pair is corresponding to the clockwise motion and the other to the 

ant i- clockwise one. 

3. Discussion 

3.1. Numerical test 

The above formalism is discussed in an idealized world. Numerical tests are done below 
to see whether the analytic result works for practical cases. First, Eqs. (13) and (14) assume 
that one can integrate observed quantities. In practice, however, observations are discrete, for 
which the integration becomes a summation. The integration and the summation could agree 
in the limit that the number of observations iV approaches the infinity. In a real world, iV 
is a large number but much smaller than the infinity. Does the above formalism still give a 
reliable answer? For investigating this point, we perform numerical simulations. According to 
the simulation for observations with equal time interval, the above formalism recovers perfectly 
the orbital parameters for N = 100. 

Next, we consider observation noises. The above formalism assumes that the observed 
points are located on an apparent ellipse. However, position measurements are inevitably 
associated with observational errors. Therefore, we perform numerical simulations by adding 
Gaussian errors into the position measurements as x — > x + Ax and y — >■ y + Ay, where Ax 
and Ay obey Gaussian distributions with the standard deviation a. We consider two cases: 
a — 0.1 (smaller case) and 0.5 (larger one) in the units of = 1. Table 1 is a list of the 
orbital parameters that are recovered by using the above formalism. See Fig. 5 for simulated 
points in a iV= 100 simulation for a = 0.1 and 0.5. For small observation error cases, the 
orbital parameters are well recovered. For large observation errors comparable to a half of 
the semimajor axis, however, the recovered angles i, 00 and ft are far from the true ones. On 
the other hand, the eccentricity ex and the semimajor axis ax are recovered better than the 
angles. The semimajor axis is overestimated, because the simulated second moments apparently 
become larger than the true ones owing to such a large dispersion. We numerically study also 
different parameter values. They lead to similar results, for which numerical tables are omitted 
for saving the space. 

Our numerical tests for the discreteness of observations and for the observation noises 
suggest that the above formalism derived for the idealized system could work in practice if 
observation noises are not so large. 

The previous analytical method cannot guess orbit parameters for such large measure- 
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ment errors a = 0.5, mostly because cosi apparently becomes larger than the unity (namely, 
the data points fit better with open orbits) (Asada et al. 2007 ). 

Up to this point, we have assumed that the orbital period is known. What happens for 
unknown binaries? We have made Fourier analyses for numerically simulated points with time 
in Fig. 5. In the time domain, the Fourier spectrum has two peaks. One peak corresponds to 
the orbital period and the other is around the artificial time step in the numerical simulations. 
This suggests that the moment approach can be applied also to unknown binary systems, if a 
Fourier analysis is adequately used to know the orbital period. Namely, the present method 
could be used to search new binary systems. 

Let us consider Cyg-Xl, a ~ 0.2 AU at 2 kpc from us. The expected angular accuracy 
in JASMINE is ~ 10 microarcsec, so that the semimajor axis of Cyg-Xl can become a direct 
observable. Other known X-ray binaries seem too faint to be observed by JASMINE. 

A Sun-like star at 20 pc with Jupiter-like planet at 1 AU could produce a wobble of 
0.001 AU, corresponding to 50 microarcsec, which must be an interesting target. 

3.2. Proper motion of the binary 

In the main part of this paper, we ignore the proper motion of the binary system in 

our galaxy. This is mostly because, for a close binary, the proper motion of the binary in our 

galaxy causes a larger cumulative displacement than the orbital motion of the component stars, 

though the orbital velocity may be larger than the proper motion. In advance, therefore, we 

know the proper motion before determining the orbital elements. For instance, it can be done 

by a comparison between the Hipparcos data and the future space astrometry. If one wishes to 

determine the proper motion (v x ,v y ) in the present formalism, however, the apparent positions 

should be replaced as x -^x + v x t and y -^-y + v y t. By averaging the observed position according 

to Eq. (13), we obtain 

1 / T 

< x > = -T obs v x + const. + O ( -=^a K 

£ V J obs 

1 / T - \ 

<y> = -T obs v y + const. + O ( ) , (57) 
^ \J-obs y 

where the terms of C^or-T^T^]) in the right hand side come from the Kepler motion. Hence 
the terms are nothing but a decaying term such as T~ b \ cos u Q b s for the eccentric anomaly u Q b s 
corresponding to t = T obs . They become negligible as T ohs T^ 1 — >• oo. Furthermore, the cumulative 
translation by the proper motion exceeds the oscillatory displacement by the Kepler motion 
for a close binary. That is, aK{T b s v)~ l <C 1, where v = J^+^- Owing to this effect, the 
decaying part due to the Kepler motion goes away for a long observation period T Q b s , say > 1 
yr. 

The parts growing linearly in the observation time T obs give the information about 
{v x ,v y ), provided T Q b s is taken as a variable in the data analysis. If T a b s 3> Tk, the above 
extraction of the linearly growing part will be possible, especially for a short-separation binary. 
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3.3. Mildly relativistic compact binary 

Finally, we mention a mildly relativistic compact binary, in which the relativistic advance 
of the periastron occurs. Rigorously speaking, we have to take account of the general relativistic 
equations. In practice, however, the above method may be applied. For instance, let us imagine 
two-year observations. For the first-year data, the periastron direction vo\ is derived. The 
second- year data tells ro 2 . The difference between the first-year and second-year directions of 
the periastron suggests the periastron advance as [w 2 — ^ 7 i)y r ~ 1 - 

For instance, the Hulse- Taylor binary pulsar shows the large periastron shift rate w = 4 
deg. per year, though the angular orbital radius is a few microarcsec. and it is beyond the 
current measurement capability. 

4. Conclusion 

This paper proposed a moment approach for orbit determinations of a close binary 
system with a short orbital period from astrometric observations alone. As an exact solution 
to the coupled equations, the orbital elements are written in terms of the second and third 
moments of the projected position that is measured by astrometry. 

The moment formalism does not replace the standard method using Kepler equation. It 
is safer to say that the present formalism is a supplementary tool for giving a rough parameter 
estimation, which can be used as a trial value for full numerical data fittings. It is interesting 
to make numerical tests of the present method. It is left as a future work. 

In the moment approach, the temporal information is smeared. Therefore, the orbital 
period cannot be determined by this method. Another method such as the Fourier analysis of 
position data with time (in Fig. 5 for example) would give a characteristic frequency that is 
the inverse of the orbital period. Fourier analyses recover the orbital period from numerically 
simulated data for Fig. 5. This suggests that the moment approach can be applied also to 
unknown binary systems, if a Fourier analysis is adequately used to know the orbital period. 
Namely, the method could be used to search new binary systems. 

We would like to thank Professor N. Gouda for useful information on astrometry mis- 
sions. We wish to thank the JASMINE science WG member for stimulating conversations. 
We would be grateful to Y. Sendouda, R. Takahashi and K. Izumi for helpful conversations on 
numerical simulations. This work was supported in part (H.A.) by a Japanese Grant-in- Aid 
for Scientific Research from the Ministry of Education, No. 21540252 (Kiban-C) and in part 
(K.Y.) by JSPS research fellowship for young scientists. 
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Fig. 1. Comparison between the second moments M xx and M yy . The orbit denoted by the closed solid 
curve has a larger variance along the x-axis, where M xx is larger than M yy . On the other hand, for the 
orbit denoted by the closed dashed curve, y components of the position have a larger scatter, where M yy 
is larger than M xx . 
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Fig. 2. Dashed curves denote two orbits: One is the semimajor axis along the x-axis. The other is along 
the y-axis. The orbit denoted by the solid curve is not distinguished by using the second moments M xx 
and M vv . The moment M xv is thus needed. 



14 



Periastron 




Actual 

Keplerian 

Orbit 



Fig. 3. Actual Keplerian orbit and apparent ellipse in three-dimensional space. We denote the inclination 
angle as i, the argument of periastron as uj and the longitude of ascending node as Vt. These angles relate 
two coordinates (x',y') and (x,y), both of which choose the origin as the common center of mass. Here, 
the x' axis is taken to lie along the semimajor axis of the apparent ellipse, while the x-axis is along the 
direction of the ascending node. 
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Fig. 4. Flowchart of parameter determinations in the moment formalism. The starting point is evaluating 
the moments from astromctric observations. 
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Fig. 5. One hundred observed points of the same source star in the x — y plane for a simulation with 
constant time interval. The parameters arc uj = 30 [deg.], ft — 30 [deg.], i = 30 [dcg.], e#- = 0.1, ax = 1.0 
and N = 100. The black curve and the gray (red in color) one denote orbits for the true parameter and for 
the mean value of the recovered parameters, respectively. Top: Gaussian errors a = 0.1. Bottom: a = 0.5. 
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Table 1. Reconstructing the parameters for numerical simulations for three different eccentricity cases as en = 0.1,0.3,0.5. 
In the table, the row cr = indicates true orbital parameters, whereas the rows a = 0.1 and 0.5 provide the recovered values 
for adding Gaussian errors (0.1 or 0.5 in the units of the true semimajor axis, respectively). For each parameter set, 100 
runs are done and the mean and the standard deviation are also evaluated. 



a 






i [deg.] 


lu [deg.] 


n [deg.] 





0.1 


1.0 


30.0 


30.0 


30.0 


0.1 


0.09454 ±0.02518 


1.011 ±0.01423 


30.02 ±2.039 


22.37 ±17.45 


29.81 ±4.205 


0.5 


0.1445 ±0.09343 


1.257 ± 0.07624 


28.87 ±7.731 


43.87 ±25.79 


40.75 ±22.62 





0.3 


1.0 


30.0 


30.0 


30.0 


0.1 


0.2936 ±0.02269 


1.011 ±0.01526 


30.04 ±2.228 


28.73 ±8.557 


30.20 ±4.687 


0.5 


0.2327± 0.1211 


1.266± 0.08215 


29.63 ±7.189 


41.26 ±25.13 


35.45 ±18.14 





0.5 


1.0 


30.0 


30.0 


30.0 


0.1 


0.4793 ±0.03103 


1.005 ±0.02253 


30.34 ±2.477 


27.53 ±8.234 


31.72 ±5.709 


0.5 


0.3065 ±0.1223 


1.233 ±0.07314 


31.33 ±5.519 


39.78 ±26.46 


40.78 ±16.08 
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